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Abstract 



An essentially exact solution of the infinite dimensional Hubbard model is 
made possible by using a self-consistent mapping of the Hubbard model in 
this limit to an effective single impurity Anderson model. Solving the latter 
with quantum Monte Carlo procedures enables us to obtain exact results 
for the one and two-particle properties of the infinite dimensional Hubbard 
model. In particular we find antiferromagnetism and a pseudogap in the 
single-particle density of states for sufficiently large values of the intrasite 
Coulomb interaction at half filling. Both the antiferromagnetic phase and 
the insulating phase above the Neel temperature are found to be quickly 
suppressed on doping. The latter is replaced by a heavy electron metal with 
a quasiparticle mass strongly dependent on doping as soon as n < 1. At 
half filling the antiferromagnetic phase boundary agrees surprisingly well in 
shape and order of magnitude with results for the three dimensional Hubbard 
model. 



PACS numbers: 71.10+x, 75.10.Jm, 75.10.Lp, 75.30.Kz. 



1 Introduction 



The Hubbard model [0J of strongly correlated electron systems has been an 
enduring problem in condensed matter physics. It is believed to at least 
qualitatively describe some of the properties of transition metal oxides, and 
possibly high temperature superconductors 0. Using standard notation, the 
Hubbard Hamiltonian reads [0] 

n= -t £ (c, n c, t , ■ (%c ln ) (1) 

<ij>,a 

+ E [ e Kt + n i,l) + U K T - 1/2) (n u - 1/2)] 

i 

where C ii<7 (C/ jCr ) destroys (creates) an electron of spin a on site % of a hyper- 
cubic lattice of dimension d, and n it<7 = C} a C iia . Despite the simplicity of the 
model, no exact solutions exist except in one dimension, where the knowl- 
edge is in fact rather complete |§. The unusual properties in one dimension 
also gave rise to discussions of whether the behaviour in 2D could be related 
to the ID case The extent that this very special limit can serve as a 
reference point for any finite dimensional model |5], || remains controversial, 
and an examination of the model from a different point of view is clearly 
needed. 

Recently, a new approach [|7|, [8], |J based on an expansion in 1/d about 
the point d = oo has been proposed to study such strongly correlated lattice 
models. In this limit the requirement of a finite total energy per site makes 
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it necessary to rescale nonlocal interactions by an appropriate power of d~ x 
[[?], |J. While e.g. spin-exchange will essentially reduce to the corresponding 
mean field theory results || , interactions like the screened Coulomb repulsion 
in the model (1) remain nontrivial even in this limit. 



In a previous short paper [|l(J , one of us presented the first method of exact 
solution of the Hubbard model in the infinite dimensional limit. In this paper, 
we provide a greatly expanded discussion of the method and the physics. 
The model in the infinite dimensional limit retains the physics expected in 
the low dimensional model, e.g. antiferromagnetism and the formation of a 
correlation induced Mott-Hubbard pseudogap in the single-particle density 
of states. In addition, we find some very interesting structures in the single- 
particle density of states (DOS), which we believe will persist in d < oo 
and lead to well defined and observable effects, which may serve as a guide 
whether the model is indeed capable of describing the essentials of transition 
metal physics. 



2 Theory and Numerical Considerations 

The limit as d — > oo is taken such that At 2 d = t* 2 = constant and it is 
convenient to choose t* = 1 as the energy scale for the remainder of this 
paper. This rescaling of the kinetic energy automatically leads to a finite net 
effective magnetic exchange dJ ~ t* 2 /U. Generally this limit provides two 
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large simplifications. First, the free DOS for near-neighbor transfer along the 
coordinate axis can be shown to acquire a Gaussian form, i.e. J?], || 

p (t) = exp(-i 2 )/v^ . (2) 

Second, the problem reduces essentially to a local one since nonlocal dy- 
namical interactions become negligible in this limit For example, the 
irreducible self energy and the irreducible vertex functions purely are local. 

This may be seen from a diagrammatic argument [|Il|] . Consider the first 
few diagrams of the single-particle self-energy for this problem as shown in 
Fig. 1. This is a real-space representation, so each electron propagator Gij 



a 




Figure 1: First few diagrams for the lattice self energy. Here, the solid 
lines represent the undressed (U = 0) electron propagators G®j(iu> n ), and 
the dashed lines represent the local Coulomb interaction U. 

scales as ~ £~l R «~ R 'j]. Thus the second order term in Fig. 1 scales as t- 3 l R ^ R i 
and, after summing over the contribution of the nearest neighbor shell, gives a 
contribution of the order d-d~ 3 ^ 2 (since 4t 2 d = 1). This contribution vanishes 
in the limit as d — > oo. A similar argument may be applied to all terms, and 



only the site diagonal self-energy survives the limit d — > oo. Furthermore, 
since the lattice is translationally invariant Ey(za; n ) = H(iu) n )8ij independent 
of i. Thus, Dyson's equation reduces to 

Gij(iu n ) = G^(iu n ) + G° ik {i(jU n )Y,{iu n )G k j{iijU n ) (3) 

k 

with the diagrammatic equation for in Fig. 1. A completely equivalent 



argument may be used for the magnetic susceptibility [12], which reads 



+ \ Xiifan, ^„)r(i^ n , iu p ; iu n )xij{iujp, iu m ; iv n ) 



Xij(iu n ; iv„) =^ e 4q - (Rl Rj) G k (^n)G' k+q (^ n + iu n ) 

k,q 

(4) 

where (3 = (/cgX 1 ) -1 , cj„ (i/ n ) are the Fermi- (Bose-) Matsubara frequencies, 
Crk(z) is the one particle Green's function obtained from (3) and T(z, z') the 
irreducible vertex. 

Usually, quantities like the self-energy and z') have to be evalu- 
ated using a direct perturbation expansion as outlined in Fig. 1. With respect 
to U this has been done for the present model quite extensively by means of 
conserving approximations including particle-hole and particle-particle lad- 
der summations fI5|. Alternatively, one of us recently proposed a perturba- 



tion theory using the transfer term t as expansion parameter |H |. However, 
as one knows from similar studies for heavy fermion systems, the validity 
and radius of convergence of these expansions may depend critically on the 
particular choice of parameters UTEj . In the following we therefore want to 
present an approach which maps the whole problem to the solution a single 
impurity Anderson model, where several reliable ways to obtain the solution 
are known. That this connection exists was realized quite early [[UJ, but 
until now it only served to set up equations for the thermodynamic potential 



ii- 

To do this let us return to equation (3). For i = j, the diagrams con- 
tributing to the local Green's function Q = Ga may be rearranged to contain 
explicitly only local processes as illustrated in Fig. 2. The undressed Green's 
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+ 



+ 
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+ 



Figure 2: First few diagrams for the local Green's function Q (double solid 
line). The diagrams have been rearranged such that all processes occurring 
on other sites than the one considered were included into an effective host 
Green's function Q represented here as a solid line. The Coulomb interacting 
is again visualized as a dashed line. 
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function represented by the solid line in Fig. 2 must incorporate the missing 
nonlocal processes and thus is the solution to a modified lattice problem 

G(iu n ) = Ga(iu n ) = Gl(iuj n ) + G° ik (i^n)^{i^n)G ki (iuj n ) (5) 

k 

where the prime at the sum indicates that k ^ i and G®Az) is the unperturbed 
propagator. This summation restriction is necessary to avoid over-counting 
of local diagrams. In reciprocal space, equation (5) reads 

G k (iu n ) = G k (iu n ) + G^{iuj n )H{iu n ) (G k (iuj n ) - G{iuJ n )) . (6) 



This may be solved for G k (z) and summed on k to yield 

1 1 



+ E(z) . (7) 



Q{z) G(z) 

The problem closes by noting two things. First, Fig. 2 is nothing but the 
perturbation expansion of an Anderson impurity model with a host Green's 
function specified by Q. In the spirit of this analogy, one may then define an 



effective hybridization strength A(z) via [18] 



A(*) = J--* + e. (8) 

Secondly, S(z) and Gu(z) are as usual related by 

Gu(z) = I dujp (uj) — - . (9) 

J z — uj — e — T,{z) 

However, Gu(z) is by definition also the propagator of the effective Anderson 
impurity problem defined by equation (8). Let us emphasize, that, although 
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we used the formal perturbation expansion of the model (1) in terms of U 
in our discussion, it does not enter the set of equations (5)-(9) at any stage. 
This means that we did not refer to any special method to solve the impurity 
Anderson model and that equations (5)- (9) are exact. 

In the present paper we should like to concentrate on the quantum Monte 



Carlo (QMC) algorithm of Hirsch and Fye [19], |2(J as a method to solve the 
impurity Anderson model, although equivalent calculations were done using 
a self-consistent perturbation scheme (NCA) for this problem [^TJ. These 
calculations generally use considerably less computational time and, away 
from some critical region at half filling, the results found were in good agree- 
ment with the QMC |22|]. However, with respect to two-particle properties 
the QMC is simply the more adequate approach, since the NCA requires the 
solution of difficult Bethe-Salpeter type equations as soon as finite Coulomb 
repulsion is considered. 

In the QMC the problem is cast into a discrete path formalism in imag- 
inary time, Ti, where T; = I At, At — /3/L, and L is the number of times 
slices. The values of L used ranged from 40 to 160, with the largest values of 
L reserved for the largest values of /3 since the time required by the algorithm 
increases like L 3 . No sign problem was encountered in the QMC process, ex- 
cept for extremely large values of U away from half filling. Even here, the 
sign problem was mild, and easily handled with the standard methods. The 
self-consistency process described above is easily employed in the QMC pro- 
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cess. Choosing an intial E, the output of the process is Q, which may then 
be inverted to yield a new estimate for E(ia> n ) and so on until Q = Gu within 
the numerical precision. Usually 4-8 iterations are required for convergence. 
Other quantities such as the unscreened local moment /i 2 =< (n-f — nj 2 >, 
the two-particle Green's functions equation (4) etc. are calculated on the last 
iteration, once convergence is reached. Typical systematic (due to finite At) 
and statistical errors were quite small, usually less than 1%. Thus, unless 
explicitely displayed, the error bars on the static quantities may be assumed 
to be smaller than the plotting symbols. 



3 Two-Particle Properties 

As already mentioned, a variety of two-particle properties may be calculated 
with this procedure. Here we want to concentrate on the antiferromagnetic 
static susceptibility, given by (4) for iv — > at q = Q = (n, 7r, • • •). Except 
for the one particle propagators discussed in the next section, the only un- 
known quantity there is the local vertex T. However, because it involves only 
local processes, it can be determined from the susceptibility of the impurity 
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Anderson model, which may be obtained from 

Xii{iLO n ,iu m ;iv n ) = XiSu n ;iv n )5 nm 

+ 7? H Xti( iuJ n, iv n )r(iuj n , iuj p ] iv n )xu(iuj p , iu m ; iv n ) 
P p 

(10) 

Here Xu is the opposite-spin two-particle Green's function matrix, 

1 rP rP rP rd . , . . , . 

Xu (iu n ,ico m ;iu n )=--55 / dri / dr 2 / dr 3 / dT 4 e- w "' (ri " T ») e -"»fa-^) 

JO JO JO JO 



(11) 

If we substitute (10) into (4), we obtain for iv — > the matrix relation 

Xq Xq Xii Xii (1^) 

where the static susceptibility is obtained by summing 

Xq( T ) = \ Z! Xq(«w„, zo; m ) (13) 

It is expected that the Hubbard model on a cubic lattice with nearest- 
neighbour transfer will exhibit antiferromagnetism at half filling for arbitrary 
values of U due to perfect nesting This transition is signaled by the 

divergence of the antiferromagnetic susceptibility xaf calculated using the 
methods described above. Results from this approach are shown in Fig. 3 
for U = 1.5 and e = 0. The logarithmic scaling behaviour is shown in the 
inset. Near the Need temperature Tjv the data fit a form xaf oc |T — T/v| y 
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Figure 3: Antiferromagnetic susceptibility Xaf(T), versus temperature T 
when U = 1.5 and e = 0.0. The logarithmic scaling behaviour is shown in 
the inset. The data close to the transition fit the form xaf oc |T — T^\ v with 
T/v = 0.0866 ± 0.0003 and v = —0.99 ± 0.05 consistent with the mean-field 
behaviour expected for d = oo . 

with T N = 0.0866 ± 0.0003 and v = -0.99 ± 0.05. This scaling behaviour 
is consistent with that of a Heisenberg model on a lattice with an infinite 
number of nearest neighbors, for which one expects the Curie- Weiss mean- 
field form for Xaf- 

The antiferromagnetic transition temperature T/v for the half-filled model 
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obtained within the current approach is plotted as a function of U in Fig. 4a. 
1.0 



0.9 
0.8 
• 0.7 
0.6 

0.15 
0.10 
0.05 
0.00^ 



(b) □ 

□ 

□ 



□ 



□ 



□ □ 



□ 



□ 



□ 



0.5 <n < 1.0 



C3- 



(a) 



O 



o 



IE 



o 



O d=infinity 
od=3 



O 



D 



O 



2 3 
U 



Figure 4: (a) AntiferromagneticT^ and (b) /i 2 =< (n-f— nj.) 2 > (calculated at 
T = Tjy) as functions ofU for the half- filled model (e = 0). In (a) the circles 
represent the values obtained from the present approach while the diamonds 
are data for a three dimensional Hubbard model extracted from ref. 24- 

For small values of U, where the local spin moment is also small, we find that 
T/v is exponentially small, consistent with perturbation theory |24j]. For very 



large values of U, where the spin moment has saturated to its maximum 
value, one expects that the transition temperature will fall monotonically 
with increasing U, T/v ~ 1/U [p^j . This is because the antiferromagnetic 
exchange dJ ~ t* 2 /U also decreases with increasing U. Thus, one expects 
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a peak in Tn(U) for some intermediate value of U as seen in Fig. 4a. For 
comparison we included in Fig. 4a also T N {U) for the d = 3 Hubbard model 



as calculated by Scalettar et al. pq| . The shape and order of magnitude 
compare very well, although the d = 3 data always have slightly larger values, 
which we think is partially a consequence of the different analytic structures 
of the free DOS in d = 3 and d = oo [26]. In Fig. 4b the unscreened 
squared magnetic moment ji 2 =< (n^ — n±) 2 >, calculated at the transition 
T = T/v, is plotted versus U when e = 0. For the half-filled model /i 2 ranges 
from /i 2 = 0.5 in the uncorrelated limit (U = 0), to fi 2 = 1 in the strongly 
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Figure 5: Critical values of U, where xaf diverges, versus filling for several 
temperatures. The result is symmetric around < n >= 1. 
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correlated limit (U — > oo). Note that the peak in T^{U) occurs near the point 
where /z 2 begins to saturate to one. Away from half filling, the divergence 
of the antiferromagnetic susceptibility is quickly suppressed. This behaviour 
is shown in Fig. 5 where the critical value of U is plotted versus filling for 
several values of (3. 

In contrast to the antiferromagnetic susceptibility, the ferromagnetic sus- 
ceptibility (q = (0,0, ■••)) is rather featureless both as a function of tem- 
perature and filling. This behaviour is shown in Fig. 6 where xf is plotted 




0.5 0.6 0.7 0.8 

<n> 



Figure 6: The Xf versus filling, when U = 4. xf was not calculated at half 
filling when f3 = 14.4 and f3 = 28.8 since xaf has already diverged here. 

versus filling for several temperatures when U — 4. For values of U as large 



13 



as eight, the ferromagnetic susceptibility never diverged. Special attention 
was placed near the region just off half filling for large U where one might 
expect Nagaoka behaviour. In this region, for U = 8, xf did show a mild 
peak as a function of filling. However, this data was prone to excessive error 
bars, and in all such cases xaf ^ Xf- 

In addition to ferromagnetic and antiferromagnetic behaviour, it is pos- 
sible that the infinite dimensional Hubbard model may exhibit incommensu- 
rate order away from half filling. Such calculations were not included here 
but are left for future study. 

By a method very similar to that used to calculate the magnetic suscepti- 
bilities, one may also calculate the superconducting pair-field susceptibilities. 
Since singlet superconducting order paramters with the same symmetry (eg. 
s-wave and extended s-wave) mix, the possible divergencies of all pair field 
susceptibilites with the same symmetry must coincide. Thus, it is only nec- 
essary to calculate one susceptibility of each symmetry. Furthermore, in the 
infinite dimensional limit, a pair-field susceptibility which corresponds to a 
symmetry orthogonal to the lattice symmetry (extended s-wave symmetry) 
has no vertex corrections. Consider a superconducting channel which has a 
form factor t k such that 



The first two diagrams in the corresponding pair-field susceptibility are shown 




(14) 



k 
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in figure 7. When d = oo the irreducible vertex has no k- dependence so the 




Figure 7: First two diagrams for the superconducting pair-field susceptibility. 
The form factor t^ determines the symmetry of the pairing. If the symmetry 
oftk is orthogonal to s-wave, then the vertex corrections vanish. 

sum over k may be performed independently on the right and left-hand sides 
of the vertex in the second term. Thus, when t k is consistent with Eq. 14, 
this term vanishes. The only diagram which remains is the simple bubble 
which can only diverge in the zero-temperature limit. This argument may 
be extended to any superconducting order parameter which is orthogonal 
to s-wave, leaving only the possiblity of superconductivity with s-wave sym- 
metry. However, for all parameters dicussed here we found that the s-wave 
pair-susceptibility remained finite at a value below the non-interacting limit. 
Thus, we conclude that there is no singlet superconductivity in the infinite 
dimensional Hubbard model. 
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4 Single- Particle Properties 

In this section we will discuss the single-particle properties of the Hubbard 
model (1) as d = oo. The most important quantity in this connection is the 
single particle density of states defined by A(u) = — 1 / 7iImGu(u +i0 + ) . Since 
the QMC method supplies us only with the imaginary time Green's function 



Gu(t), we used the maximum entropy procedure |27], |2£| for providing 
its analytic continuation. In most cases results from a self-consistent per- 
turbation theory for the Anderson impurity (NCA) ||2l| , [30fl were used as a 
default model for this procedure. This is a nearly ideal combination of meth- 
ods, since the NCA becomes exact at high frequencies where the QMC data 
contains little information, and the analytically continued QMC data is es- 
sentially exact at low frequencies. A detailed comparison of NCA and QMC 



results will be presented elsewhere f22fl . When a pseudogap was present in 



the one-particle excitation spectrum, we found that the NCA equations failed 
to converge due to numerical instabilities. In this limit second-order pertur- 
bation theory in U was used to provide a default model for the maximum 
entropy process. 

With either default model, it was extremely important to carefully con- 
trol both the systematic and especially the statistical errors of the Monte 
Carlo process. Thus, before analytic continuation, the self consistent process 
discussed above was allowed to fully converge. Then several runs were per- 
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formed to obtain a highly accurate average value of E. The self-consistency 
was then turned off, and several runs using this value of £ were performed 
to calculate a highly accurate value of Q (r) which was then analytically con- 
tinued. 

Results for the half-filled model are shown in Fig. 8. Here A(u) is plotted 
for several values of U when (5 — 7.2, e = 0, and < n >— 1. As the Coulomb 
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Figure 8: Evolution of A{uj) in the paramagnetic phase of the half-filled model. 

repulsion U is increased from zero the broad central peak in the spectrum 
becomes narrower and is gradually suppressed while at the same time two 
side bands build up which are carrying the majority of the spectral weight. 
This kind of behaviour has also been found in different kinds of perturbation 
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theory [JTJ. As discussed below, the central peak can be identified with 
a Kondo-like resonance, and the side peaks as usual with charge transfer on 
and off the site. As U continues to rise, the spectrum begins to develop 
a pseudogap at zero frequency when U > 3.4 « Uc- The "critical" value 
Uc increases with increasing temperature, and for U > Uc the pseudogap 
grows linearly in U. This structure is identified as a pseudogap, rather than 
a gap, since the zero frequency density of states was never identically zero. 
Rather at ou — it was about 4 orders of magnitude smaller than that of the 
surrounding features. Indeed, this is consistent with an argument of Krauth 
and George who argue that in this model the zero frequency density of 
states is finite for all T, even at T = 0. 

The evolution of density of states and the filling as a function of e is shown 
in figure 9. As e is increased from zero, the filling decreases very slowly at 
first, consistent with the presence of a pseudogap. The filling begins to change 
rapidly once e is increased to roughly half the gapwidth. The change in filling 
is enhanced by the presence of a narrow resonance near zero frequency which 
replaces the pseudogap as soon as the model is doped away from half filling. 
However, there is a vestige of the Mott pseudogap at higher frequencies. 

Perhaps even more interesting is the temperature dependence of the DOS. 
In Fig. 10 the DOS for U = 4 close to half filling (n ~ 0.94) is displayed for 
some typical values of T. While the peaks related to charge excitations are 
nearly temperature independent, the resonance at the chemical potential is 
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Figure 9: Evolution of A{uj) as a function of filling when (3 = 7.2 and U — 4. 
The filling versus e is shown in the inset. 

very sensitive to the temperature and evolves roughly as ln(T) when going 
from /3 = 7.2 to /3 — 14.4. It becomes more pronounced at lower tempera- 
tures, and apparently disappears completely once the temperature is raised 
to about twice the width of this feature. As shown in the inset in Fig. 10, the 
growth of this resonance is correlated with the reduction of the screened local 
moment Txu(T). In connection with the fact that eqs. (5)- (9) represent an 
effective Anderson impurity problem, it is evident that this behaviour must be 
interpreted as a Kondo-like screening of the local moment and consequently 
the resonance at the Fermi surface can be identified with the Abrikosov-Suhl 
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Figure 10: Evolution of the Kondo-like feature with T. As shown in the in- 
set,the growth of the resonance is correlated with the reduction of the screened 
local moment T\u- 

resonance known to accompany the former. With the reasonable assumption 
E(i0 + ) -> as T -> 0, we expect A(0) -> 1/y/n @. Using as the "Kondo"- 



scale the temperature, where A(0;T K ) = \A(0, T — > 0), we find Tk ~ t*/8 
in this case. 

Since this Kondo-like feature is pronounced, and occurs right near the 
Fermi surface, one would expect it to have a strong effect upon a number of 



properties of the system like resistivity and specific heat p2[ . One simple way 
to obtain an idea how large this influence may be is to inspect the effective 
mass of the quasiparticles in the system. In Fig. 11, the mass enhancement 
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factor 



a 1 (T) = 1 — (Im£(zo;o)) /uiq 



(15) 



is plotted versus temperature when U — 4 for some typical values of e. Here, 
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Figure 11: Mass enhancement a X (T) = 1 — (ImS(zo;o)) A^o versus tempera- 
ture. 

iuiQ is the lowest Matsubara frequency, and the wave function renormaliza- 
tion factor is given by the zero temperature limit of a(T). As shown in 
the figure, the mass enhancement factor increases like ln(T) until the tem- 
perature falls well below the Kondo-scale, in which case it saturates. The 
behaviour a" l (T) ~ ln(T) has been identified @] as one of the signatures 



of marginal Fermi liquid behaviour which has been suggested as a possible 
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phenomenological picture for the two-dimensional limit of the model (1) [|35 
However, it is clear that the infinite-dimensional model close to half filling is 
not marginal, since a -1 saturates to a finite value at low temperatures, cor- 
responding to a usual Fermi liquid with an enhanced effective mass, which 
is strongly dependent upon filling. In the case when e = 1.0, the filling is 
roughly < n >^ 0.94, and a _1 (0) ~ 7.5, corresponding to a mass enhance- 
ment of m* /m « 7.5 ~ T£ . Far away from half-filling, it is about one, and 
as half filling is approached it increases dramatically. 



5 Conclusion and Interpretation 

The method described in this paper has reduced the infinite dimensional 
Hubbard model to a self-consistently embedded Anderson impurity problem. 
Thus for large U the qualitative features of the density of states have an 
obvious interpretation in terms of the Anderson model spectrum. The upper 
and lower peaks found correspond to charge transfer on and off the local site. 
In addition a small resonance at or near the chemical potential occurs which 
to our opinion must be interpreted in terms of an Abrikosov-Suhl resonance 
accompanying the screening of the local moment as observed in Txu{T). As 
well known from the physics of the Kondo effect, this formation of a bound 
state causes a significant mass enhancement for the quasiparticles in the 
system as T — > 0. 
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From the data presented here, it is possible to obtain a very good repre- 
sentation of the phase diagram of the infinite dimensional Hubbard model. 
For example, the phase diagram for half-filling is illustrated in Fig. 12. At 
0.4 
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H 0.2 

0.1 

0.0 

1 2 3 4 5 

u 

Figure 12: Phase diagram for the Hubbard model in d = oo at half filling. The 
shaded region is centered upon the location of the continuous metal-insulator 
transition. The extension of this boundary into the antiferromagnetic region 
is represented by the dashed line. 

temperatures above T N we find a paramagnetic Fermi liquid state with more 
or less enhanced quasiparticle masses. This evolves gradually into a Mott- 
Hubbard phase which is stable beyond some critical value of U. This "insula- 
tor" phase is identified with the occurance of an exponentially small density 
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of states at zero frequency and is quickly suppressed upon increase of T, 
because it can only be stabilised when the charge excitation peaks are well 
separated on the scale of variation of Fermi's function. Since the zero fre- 
quency density of states is always finite, no real phase transition (i.e. one 
with an identifiable order parameter) from the metal to the insulator has 
occurred. However, when the zero-frequency density of states becomes expo- 
nentially small, the dynamic and and transport properties [^] of the system 



change radically in a manner consistent with the change from a metal to an 
insulator. At low temperatures the system is always antiferromagnetically 
ordered with a symmetry-induced gap in the one-particle DOS. Upon doping 
(not shown), the Mott insulator phase vanishes, and the AF region shrinks. 

If we artificially extend the metal-insulator boundary into the antifer- 
romagnetic region (the dashed line in Fig. 12), then we find an interesting 
reentrance: the critical U c increases with decreasing temperature. This reen- 
trance is not relevant to the present unfrustrated model due to the antiferro- 
magnetism; however, it is relevant to models with sufficient lattice frustration 
so that the antiferromagnetic transition is suppressed. Such is the case if we 
include a hopping term to the second next neighbor along each of the co- 
ordinate axes equal to the nearest neighbor hopping t. In fact, with this 
modification, the unperturbed density of states in the d = oo limit remains 



Gaussian, with a renormalized energy scale [33]. Our present calculation is 



relevant to such a model. In our opinion this reentrance must be attributed 

24 



to a competition between the Kondo-effect, which never completely vanishes 
due to the finite DOS at the chemical potential, and the Mott-Hubbard 
phase. Actually, from inspection of the ground state energy we find that the 
Abrikosov-Suhl resonance is always favoured when T — > and U large. We 
emphasize that this behaviour is strongly related to the fact that for d = oo 
one always must have ^4(0) > 0, which leads to a small Kondo-effect that 
eventually kills the Mott-Hubbard transition as T — > 0. Without symmetry 
breaking due to a magnetic phase transition we expect a similar situation in 
finite dimensions and T > fll4|l . However, a smaller value of U c as T —>■ 
must be anticipated in this case, because ^4(0) should decrease exponentially 
with temperature for sufficient large U. 

In conclusion, we have presented a set of self-consistent equations for 
the Hubbard model which allows one to calculate the properties of strongly 
correlated systems in the limit of infinite dimensions. Using a quantum 
Monte Carlo procedure to solve these equations we have shown that the 
model displays the expected antiferromagnetism when half-filled, and that 
the single-particle density of states displays a correlation pseudogap. The 
importance of this result is that it allows an essentially exact solution of 
the d = oo Hubbard model in the thermodynamic limit. Thus, solutions 
now exist for the model in two limits, d = 1 and d = oo. With respect to 
the various phenomenological models it is surely noteworthy that the feature 
found at the Fermi level provides to some extent a natural way to obtain 
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unusual properties. Finally, while this method is discussed in the context of 
the Hubbard model, it could be applied to other strongly correlated models 
(i.e. the periodic Anderson model) just by changing the form of the undressed 
Green's functions in (5)- (9). Work on this problem is in progress. 
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